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ABSTRACT 


Results from continued research on selected parameters to minimize community 
annoyance from airport noise are reported. First, a review of the initial 
work on this problem is presented. Then the research focus is expanded by 
considering multiobjective optimization approaches for this problem. A 
multiobjective optimization algorithm review from the open literature is 
presented. This is followed by the raultiobjective mathematical formulation 
for the problem of interest. A discussion of the appropriate solution al- 
gorithm for the multiobjective formulation is conducted. 

Alternate formulations and associated solution algorithms are discussed 
and evaluated for this airport noise problem. Selected solution algorithms 
that have been implemented are then used to produce computational results 
for example airports. These computations involved finding the optimal operating 
scenario for a moderate size airport and a series of sensitivity analyses 
for a smaller examp ’’e airport. 


CONTINUED RESEARCH ON SEV 'ITED RARAMETERS TO MINIMISE 


COMMUNITY ANNOYANC’ FROM AIRPORT NOISE 


1. INTRODUCTION 

This report describes the continued research on selected control param- 
eters for minimizing coimnunity annoyance with airport noise. This research 
provides a means of assessing the relative importance of problem parameters 
and their impact i:i terms of annoyance on conmiunities subjected to airport 
noise. 

The second section of this report briefly reviovi?s the work, performed 
on this problem during the initial one year's investigation. Section three 
provides a literature review for Multiobjective Optimization techniques. 
Alternate solution algorithms and formulations are discussed in the fourth 
section of this report. Computational Results for a moderate size airport 
are reported upon in section five once the mathematical formulation is made 
available. The concluding section of this report presents a series of 
sensitivity analyses conducted during the course of this study. 

2. INITIAL STUDY ON SELECTED PARAMETERS TO MINIMIZE COMMUNITY ANNOYANCE 

FROM AIRPORT NOISE 

The initial research on the airport noise problem first concerns itself 
v.itL the development of mathematical models that functionally represent the 
problem of interest. These models may be classified as corresponding to 
nonlinear integer mathematical programming problems. The objective function 
in every case represents the minimization of a noise annoyance measure repre- 
sented as a function of airport activity. The controls (decision variables) 
manipulated to achieve this objective were aircraft fleet mix composition 
for the airport along with runway assignment as a function of the time of 

- 1 - 


d 


day. ConBiderad aa constraints were- 


1. The demand for f light services 

2. The fleet mix cotiipnsition flexibility 

3. A limitation on daily flights 

i\. Restrictions on the use of specified runways as a 
function of type of aircraft and/or time of day 
The constraints as formulated were linear in terms of the decision 
variables. 

The solution algorithms investigated as methods to obtain solutions for 
the mathematical models developed were separable programming, relaxation 
and partitioning procedures. With relaxation and partitioning one solves 
special structured problems, for which efficient solution procedures exist, 
in an iterative fashion. Separable programming Involves successive linear 
approximation to obtain solutions to the nonlinear problem of interest. 

These solution algorithms were applied upon several applications and optimal 
(i.e., minimum annoyance) operating scenarios were obtained. Xt was deter- 
mined that the relaxation technique that involved the solution of succes- 
sive linear programs was the most efficient algorithm for the type of air- 
port problems faced. 

3. MULTIPLE OBJECTIVE I'UNCTIONS 

It is recognized that working singularily towards the objective of 
minimizing annoyance from airport noise can degrade other potential goals 
such as fuel consumption and aircraft time delay. Because of this the 
additional objectives of fuel consumption and time delay of aircraft are 
now considered. 




\ 


^ 4t :A rf*. 


I 






I 

I 

I 

! 

I 

i 


I 


I 

\ 


3.1 MULTIPLE OBJECTIVE FUNCTION FORl'rLATION 

There are various types of air l••rrainal time delay. Reference (11) 
identifies four types of such delays. Namely they were identified as 
holding, path-stretching, early descent/speed control and procedural routing. 
The delay of interest in this research is that of procedural routing. 

Probably the most common type of procedural routing is designed to comply 
with local noise abatement policy. It has been shown in (11) that the ad- 
ditional fuel consumption due to procedural routing is a linear function of 
the associated time delay. Time delay can be determined for any aircraft 
knowing its speed and the trajectory selected as well as the minimum time 
path. Notation: 


'‘idj 

'"ikdj 

‘^idj 

^ikdj 


T 


F 


ND 

^ikdj £ 


^'idj 1 


the time delay for the i — type aircraft arriving from 
direction d and utilizing trajectory j. 

til 

the time delay for the i — type aircraft of stagelength 
k departing for direction d and utilising trajectory j. 

the additional fuel consumption for the i~ type air- 
craft arriving from direction d and using trajectory j. 

the additional fuel consumption for che i~ type aircraft 
of stagelength k departing for direction d and utilizing 
trajectory j. 

the total amount of time delay allowable. 

the maximum amount of additional fuel consumption allox^ed. 

the number of directions considered. 


number of departures of type i aircraft with stagelength 
k going in direction d utilizing trajectory j during period 

£. 


number of type i arrivals from direction d utilizing tra- 
jectory j during period £. 


The alternative objective of minimizing the total time delay may be for- 
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The objecuive of minimizing the total amount of additional fuel consumed 
may be represented by 
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The additional constraints considered were, 


0„ < T 

(3.3) 
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Constraint (3.3) restricts the amount of time delay allowed x^here as 
(3.4) specifies an upper limit on the additional fuel which may be cpnsumed. 

To complete the multiple objective formulation equations the original 
model (Appendix A) xwuld be changed to account for the extra dimension 
(direction) on the decision variables. 

A multiobjeccive solution algorithm has of yet not been implemented 
since the noise objective function xras of primary concern to the users of 
this research. An alternative to a multiobjective solution algorithm would 
be to treat the objectives of time delay and additional fuel consumption 
x^ith constraint equations (3.3) and (3.4) and then perform extensive 
sensitivity analyses by varying the right-hand-side of these equations. A 
fuel minimization computer run xjas accomplished and the results reported 


on in Appendix E. 


3.2 LITERATURE REVIEW 

The consideration of the alternative objective functions of minimiziinp, 
fuel consumption and minimixing passenger and or freight delay leads one to 
consider the potential of vector optimization solution procedures. The 
following literature review is meant to provide background material for the 
consideration of multiple objective optimization. 

Hwang (12) has classified multiple objective decision methods (MODM) 
into four categories based on the availability of preference information 
from a decision maker. They arc, 

1) no articulation of preference Information 

2) "a priori" articulation of preference Info-'mation 

3) progressive articulation of preference information 
(iterative methods) 

4) "a posteriori" articulation of preference information 
(generating techniques) 

The last three categories are of primary interest and will be discussed 
in greater detail. 

3.2.1 "A Priori" Articulation of Preference Information 

Those methods depend on the ordering or prioritizing of the objectives 
of interest. This helps to reduce the size of the solution set which must 
be considered. The primary difference between the various methods that 
belong to this MODM category is in the nature of the preference information 
required. Several classes of methods that belong to this category are now 


Utility Function Method 

This method ia based upon obtai'iing cardinal information from the 
decision maker in the form of specific objective preference levels. A 
satisfactory solution can be obtained if the proper utility function is 
available and is used in conjunction with the decision maker's preferences. 
The prime difficulties in using such methods are the specification of the 
appropriate utility function and the computational burden that increases 
exponentially with the number of objectives considered. 

Goal Programming 

Goal programming methods seek to minimize the absolute deviations of 
the objectives from predetermined targets or goals. Originally proposed 
by Charnes and Cooper (3) many different approaches have been developed. 

Lee (13), Ignizio (12) and Cohen (4) describe many such methods. 

The various approaches differ in the formulation of the surrogate 
objective function. Advantages cf Goal Programming is that it is relatively 
efficient computationally and requires less effort by the decision maker 
in specifying preference information. The major disadvantage of Goal Pro- 
gramming is that any solution determined is a function of the goals set 
by the decision maker. Unless the decision maker is careful in establish- 
ing a priority structure one can arrive at an infeasible or dominated 
solution. 

Electro Method 

The Electre Method, proposed by Roy (IS) is based upon an "out- 
ranking" relationship for non-inferior solutions. A partial ordering of 
the solutions is obtained based upon two conditions; a concordance con- 
dition and a discordance condition. The function of the concordance con- 
dition is to indicate a tolerance for errors in comparing candidate solu- 
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tions, The diacorUanco eondition is » icd to indicate which Golutions need 


to be compared. 

3.2.2 Iterative Methods 

AIJ- of these methods basically utilize two steps in obtaining solutions 
I’irst a technique is employed to generate a solution. Secondly, methods 
are used to elicit the decision makers reaction to the given solution and 
then use this reaction to modify the problem by a modified priority 
structure. Two such methods are now reviewed. 

Surrogate Worth Trade-off Method (SWT) 

This method proposed by Haines, Hall and Freedman (10) consists of 
(1) generation of candidate (nondominated) solutions which are used te 
form trade-off functions and (2) the search tor a preferred solution using 
a surrogate worth function. The trade-off function is formed from the rel- 
ative trade-offs of marginal Increases or decreases from any two objective 
functions. The surrogate worth function is the decision maker's assessment 
of how much one prefers trading one objective for another. The trade-off 
functions are generated for any two objectives assuming that all remaining 
objectives are fixed. This means it provides limited information to the 
decision maker about the nondominated solution set. 

Step Transform Evaluation Method (STEM) 

The use of Linear Programming (LP) with Multiple Objective Functions 
Ste.p Method has been proposed by Benayoun and others (1) , This method 
involves, 

1) construction of a pay-off table after optimizing 
with p objective functions separately 


2) oolving an LP w' idi io "nearoat" in the MINIMAX 
case to an idea! solution. 

3) specifying a trade-off roiatlonohip by the deeision 
maker. 

Other interactive methodo arc described by Geoffrion <9), Zionto and 
WaXliniUB (20). 

3.2.3 Generating Techniciucs 

Generating Techniciucs involve the generation of the nondoroinated solu- 
tion set. These methods seek to develop an exact or approximate repre- 
sentation of the entire solution set. Because of this they are computa- 
tionally intensive .and feasible only for problems with a small number of 
objectives . 

The Parametric method proposed by Gal and Nedoma (8) and Zadeh (19) 
and the Epsilon-Constraint method due to Marglin (14) ore examples of 
Generating Techniques that have been applied to multiple objective decision 
problems . 

Another approach sometimes used is to relegate all objectives 
except the one considered most important to constraint equation status. 

The problem is first solved with the single selected objective function. 
Post optimality analysis is then employed by varying the right-hand-sides 
of the constraints that originally corresponded to alternate objective 
functions. Thus, an entire family of solutions may be generated as a 
function of the alternate objective functions. It is this approach that 
is presently being implemented to deal with the alternate objectives of 
additional fuel minimization and passenger/freigbt delay. 


In order Co ncleec u oolucion o 'Qtegy for a given ir.ulti- 
objective problem xt io neccoonry to • omparc the effectivoncoo of potential 
oolucion techniqueo for the opecific problem of intcroot. Normally one 
would evaluate candidate solution algorithmo on the basis of criteria ouch 
as 


1) computational burden 

2) quantification of tradeoffs 

3) availability of decision making information 

Such research has been initiated and v^ill continue over the coming year. 

h. ALTERNATIVE FORilULATIONS AND SOLUTION ALGORITHMS 

The research into the design and development of efficient solution 
algorithms for minimizing conununity annoyance from airport roise has been 
continued, The reason for this is that the simplest and most efficient 
solution algorithm possible is desirable. However, such an algorithm 
should produce superior results for the problem of irvterest. 

The solution algorithm implemented for the airport noise problem 
(see Appendix B for a description) does not guarantee a global optimum. 
However it does appear capable of always generating very good solutions 
for the airports studied so far. In some iustances the use of the selected 
solution algorithm would be quite costly and time consuming. Because of 
this several alternatives to the implemented solution algorithm are pro- 
posed, discussed and limited computational experience reported. 

For airports with a large fleet of many different servicing aircraft 
that utilize many different departure and arrival trajectories the number 
of decision variables and constraint equations may become prohibatively 
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large. In ouch cocco a olmplifled p >blem may be solved. The simpli- 
fication is that the distinction bef ‘cn aircraft types is dropped and ail 
operations are for a "typical" aircraft. The corresponding mathematical 
model is of the same form as the original model (see Appendix A) except 
that the type of aircraft designation is removed . Therefore the decision 

“ number of departures utilizing trajectory j 
during time period £. 

“ number of arrivals utilizing trajectory j 
during time period £,. 

For comparison purposes an airport with the following characteristics would 
have 2510 decision variables for the original model whereas it would in- 
volve only 60 variables with the simplified model. 

14 different types of aircraft 

15 " trajectories 

2 " time periods 

5 " departure stage lengths 

The use of this simplified model would become even more attractive if it 
imperative to consider nonlinear constraints or problems with more than 
two time periods. For example if constraint (A. 3) were to be considered 
for an actual airpoi’t exercise we could not use the original solution 
algorithm since it relies on all constraints being liin^ar. Another advan- 
tage associated with this simplified model is that nonlinear optimization 
solution algorithms might be successfully applied if the number of vari- 
ables and constraints were not excessive. 

During recent years the parametric transformation technique know 
as the Augmented Lagrangian (ALAG) Penalty Function Technique has gained 
recognition as one of the most effective type of methods for solving con- 


variables are 
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trained minimization problems. In ti>j opinion of many researchers iii chis 
field, the AhAG penalty function teoiniquo is the best method available for 
solving problems with nonlinear constraints in the absence of special struc- 
( 2 ) , The disadvantage of the method are negligible and the advantages are 
strong, especially the lack of numerical difficulties and the ease of using 
the unconstrained minimization routine. The method has global convergence 
at an ultimately superlinear rate, the computational effort per minimization 
falls off rapidly, initial starting point need not be feasible and the func- 
tion is defined for all values of the parameters ( 2 ) , 

The ALAG penalty function technique is a balance between the classical 
penalty function technique and the Lagrangian primal-dual method which are 
both parametric transformation techniques. The design of this method was 
motivated by efforts to overcome the numerical instability of the penalty 
function technique near the solution (12), (18) and attempts to eliminate 
the "duality gap" in nor-’envex programming ( 6 ) . This ALAG penalty function 
nonlinear optimizing algoritlim has been Implemented to obtain solutions for 
the simplified model described previously. The documentation for this im- 
plemented solution algorithm is presented in Appendix C, Computational 
results for this algorithm are presented in Appendix D. 

Another alternative solution algorithm would be to utilize a linearized 
version of the nonlinear objective function. This would amoimt to stopping 
after step 1 of the solution algorithm reported upon in the first years 
report (see Appendix B) . The advantage of this approach is that one is 
solving Linear Programming problems for which very efficient solution tech- 
niques exist. The main disadvantage is that the surrogate objective func- 
tion is just a linearized representation of the nonlinear objective func- 
tion. Computational results, thus far, indicate that this surrogate objective 
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' function is an effective substitv.te lor the nonlinear objective function. 
That is, it doesn't require near the computational effort and y+toduces 
solutions nearly as good as the nonlinear algorithm. 

Another useful formulation of the community response objective func- 
tion would involve the goal of minimizing the maximum annoyance experienced 
by any population segment within the airport vicinity. The motivation for 
such a goal is that the annoyance due to airport noise is shared in a some- 
what uniform manner by the effected population segments. 

Several possibilities that have been considered are minimize the max- 
imum, 


1. noise impact index contribution from the various 
population segments 

2. noise exposure forecast (NEF) or day-night level 
^^dn^ for the population segments 

3. population weighted NEF or population weighted 

an 

The formulation of such goals can be accomplished by letting 

n = the number of population segments 

= response index for population segment k. 

then 


minimize {maximum ^^2' 


. , RI ]} 
n 


which may be written as 


minimize Z 

subject to; Z ^ ^ ' 1>2 n 

where Z = maximum RI 2 ’ 

An objective function utilizing a population weighted NEF has been for 
raulated ead run for the example airport number two described in ( 7 ) . The 
computational results are presented in Tables FIO and Fll of Appendix F. 


5. APPLICATION AIRPORT FOR SOLUTION PROCEDURE 

This section presents an example, medium sized airport to which 
the suggested solution procedure has been applied. The analysis contained 
herein is limited to commercial airline traffic for the particular 
airport of interest. A series of tables and figures are first presented 
defining parameter values needed in constructing the appropriate mathe- 
matical formulation. The profiles utilized are the standard ones described 
previously in (7) . The computational results obtained from application of 
the solution algorithm to this example airport are displayed in Table 5.7. 


TABLE 5.1 


Demand for Incoming Flights 


Runway 

Tracks 


NiRht 


(13)12 R 


(14)30 L 


















































TABL' 5.5 


Population of Areas in t icinity of Example Airport 


Atqix 

Pop 

Area 

Pop 

1 

20078 

34 

879 

2 

16176 

35 

2127 

3 

8022 

36 

1161 


8732 

37 

2591 

5 

11887 

38 

7447 

6 

12317 

39 

4294 

7 

4987 

40 

3585 

8 

5822 

41 

4401 

9 

7053 

42 

8060 

10 

19680 

43 

6075 

11 

9494 

44 

13940 

12 

3579 

45 

9578 

13 

5596 

46 

10334 

14 

946 

47 

6659 

15 

96 

48 

17991 

16 

8918 

49 

17268 

17 

5339 

50 

6714 

18 

9475 

51 

3048 

19 

532 

52 

8144 

20 

8874 

53 

13093 

21 

3799 

54 

5193 

22 

3996 

55 

5359 

23 

453 

56 

21192 

24 

1250 

5/ 

13785 

25 

1570 

58 

5640 

26 

4109 

59 

16827 

27 

5833 

60 

17408 

28 

6908 

61 

7977 

29 

107,16 

62 

15239 

30 

11813 

63 

10034 

31 

2399 

64 

36311 

32 

15610 

65 

10852 

33 

10661 
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TA.'^E 5.6 


Typer, of Mrci ifts Considered 


— — — "" 

Aircraft Type 

Number of Stagelengths* 

DC-9-32 

3 

B727-200 

4 

DC-8-55 

4 

L-1011 

4 


f 


I 


* The stage length of a departing aircraft is a measure of the distance 
to the. next destination. Stage lengths 1(0-500 nautical miles), 
2(500-1000 nautical miles), 3(1000-1500 nautical miles) and 4(1500- 
2500 nautical miles) 


IL. 


.■it 
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I TA '.E 5.7 

i COMPUTATIONAL RESULTS FOR A PLICATION AIRPORT (DEPARTURES) 

f 

Nil 0.17723 













6. SENSITIVITY ANALYSES 


It Bhould ue ro^ogniKed that th>‘ winimisation solution to the noise 
problem under study is optimal only with respect to the specific model being 
used to represent the real problem, and such a solution becomes a reliable 
guide for action only after it has been verified as performing well for other 
reasonable representations of the problem as well. Furthermore, the model 
parameters (particularly the right-hand-side of constraint equations) are 
set as a result of operating policy and these decisions should be reviewed 
after seeing their consequences on what can be achieved. 

For these reasons alone it is Important to perform sensitivity analyses 
to investigate the effect on the optimal solution of the various problem 
parameters taking on a range of possible values. At times there will bo 
some parameters that can bo assigned any reasonable value without affecting 
the optimality of the problem of interest. However, thero are usually other 
parameters which if perturbed in value would yield new optimal solutions. 

This is particularly serious if the change in a problem parameter results 
In a substantial change in the value of the objective function. Therefore, 
the basic goal of the on-going sensitivity analyses is to identify those 
particularly sensitive parameters, so that special care may be exercised in 
estimating them and in selecting solution algorithms that perform well for 
the most likely values of such parameters. 

The sensitivity analyses conducted as a continuation of research into 
the noise problem ares 

1. V.ir! . .. Of: in the demand for flight services 

2. Fleet mix composition changes 

3. Changes in operations by time period 

4. Variations in objective function formulation 
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Thofio oemslcivity analyeoo conoiderc • and conducted for a particular air- 
port arc reported on in Appendix fheso results are sunimaric-cd in Tables 
6.1 and 6.2. 

For eoraparison purposes tbo results from Tables E2-E3 should be grouped 
together. These Tables display the results for variations in demand for 
flight services. It is found that in a vast majority of the cases that 
northwest (NW) through easterly (K) flights are assigned track 9 
easterly throt.gh southwesterly (SW) " " " " 6 

southwest through nortlwestorly " " " " 3 

This is one indication that the assignment of aircraft to tracks is some- 
what insensitive to variations in demand lor flight services. 

Tables K6-E9 results may also be grouped together for comparison pur- 
poses. These Tables corresponded to the use of only one type of aircraft 
for all commercial service. Those cases display more variation in assign- 
ment, howe.ver, somewhat general observations can be made. 

For departures 

NW through E flights are assigned track 9 

E " SW " " " " h 

SW " NW " " " " 1 

For approaches 

NW through E flights are assigned track 9 
E " SW " " " tracks 4 and 7 

Even though greater variation is displayed in the results from Tables F,6- 
E9 (due entirely to variations in approach assignments) it is conjectured 
that the assignments for this example are relatively insensitive to the 
type of aircraft considered. The reason for this is that there are a great 
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TAB *. 6.1 


HUMMAKY OF SENSITr-nT ANALYSES RESULTS 
FOR APPROACHES 


I 



Entries in Table are percentage use 
for each period. 

i 


! 



I 
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SU^2■^ARY OF SENSITl. f.TY ANALYSES RESULTS 
FOR DEPARTURES 


FROM TABLE 


Tracks 


1 (day) 


2 (day) 




E7 E8 







'1 

21 

21 

21 

21 

13 








8 

22 

18 

20 






57 

61 

54 

56 

56 

56 

56 

50 

50 

50 


50 

SO 

SO 

50 


2.5 


2 






2.5 


9 





23 

16 

21 

22 

23 

23 

23 

23 

50 

50 

50 


50 

50 

50 

50 


Entries in Ta e arc percentage use 
for each period. 














many alternative solutions that give lust as good or nearly as good an ob- 
jective function value as the optimal ty selected solution. 

Even though the results seem to indicate that the optimal solution to 
the airport noise problem is insensitive to the variations considered, for 
the one example airport, this should not lead one to conclude this would be 
true of all airports or all possible variations. The primary reason for 
the insensifivity displayed in the reported results is that the capacity 
constraints for the approach and departure trajectories were never binding 
in any of the cases considered. For a larger busier airport these capacity 
constraints might bo ninding for some flight service scenarios and one 
would expect greater variability in sensitivity analyses. If an annoyance 
metric; were available chat accounted for human annoyance as a function of the 
time frequency of flights and this metric was utilised to obtain noise 
optimal flight schedules then one would expect greater variation in sensi- 
tivity analyses. Tables E6 through E9 provide the corresponding computa- 
tional results. 

Table ElO describes the results obtained when the objective function 
was changcid from the one originally given in equation (El) to that of mini- 
mizing the maximum (mlnimax) annoyance for the various population segments 
considered. The use of such an objective formulation should result in a 
dispersion of annoyance more uniformly to the population areas surrounding 
the airport. The solution of a population weighted NEF minimax objective 
function is given in Table Ell. 
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APPENDIX A 


AIRPORT NOISE MINIMIZATION MODEL 


Reference (7) 


A.l Notation 


CONSTANTS 


A 

i 

j 

k 

P 

R 

NA 

ND 

NI 

NJ 


NK 


NL 


NR 

NV 


N 


A 


» area dt iignation 
=> type ol aircraft designation 
= ground track designation 
» stage length designation 
= time period designation 
= total affected population 
= runway number 
«= number of areas 

a number of runways for departure aircraft 
a number of types of aircraft 
a number of different ground tracks 
a number of stage lengths 
a number of time periods 
a number of runways 

a number of runways for arriving aircraft 
a limitation on number of operations that use ground 
track d 

a limitation on number of operations in time period P 
a limitation on number of operations for runway r 
a limitation on number of operations with stage length 
a limitation on number of airport operations per day 
a population in area A 

a set of ground tracks associated with runway R 
a community response index critical value for area A 

a number of type i aircraft available for take-offs 
during time period i 



I-' 


I' 

I 

I 

I 


i 

I 

L 

I 

[, 

I 

I 


number of typ • i aircraft available for landing 
during time p Tiod i 

^iBdp number of type t aircraft with stage length s with 

ground track d that are required during time period p 


NOISE PARAMETERS 


EPNL 

EPNL 


IjA 


EPNL 


ikJA 




NEF 


W(RIa) 


° affective perceived noise level 

= effective perceived noise level for arriving aircraft 
of type i corresponding to ground track j CKperienced 
in area A 

« effective nerceived noise level for departing aircraft 
of type i with stage length k corresponding to ground 
track j experienced in area A 
“ A-Weighted sound pressure level 
“ day-night level 
“ equivalent sound level 
a noise exposure forecast 
“ community response index for area A 
“ weighting factor as a function of conmmnity response 
index for area A 


VARIABLES 

^ikj a 



« number of departures of type i aircraft with stage 
lexT,gth k utilizixig trajectory j during period a 
= number of arrivals of typo i aircraft utilizing 
trajectory j during period a 








EPNL 


NA 


NI NK NJ 


minimise [ K {P,/P)(10) lo { K Z T. (10 — 
A°1 ^ i'^l k-=^l j«l 
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(A.l) 


This is the form of the objective function used in solving the example 
problems reported on in Section IV. 

A, 2 Constraints 

Airport authorities seeking to reduce noise and still service passenger 
demand for their airports may impose various related operating constrants. 

A. 2.1 PH 'ht limitations 

For any given airport flight limitations may be imposed for many differciit 
reasons and consist of many different constraints on aircraft operations. The 
flight limitations of primary interest in this research are those operational 
constraints that may be imposed in an attempt to improve the noise enviroment 
around an existing airport. The operational flight constraints formulated in 
this research are: 

1. The annoyance as measured by a specified community 
respm^se index may not exceed critical values for 
specified communities in the airport vicinity. 

2. The number of certain types of aircraft that may 
operate into and out of a given airport cannot ex- 
ceed some upper limit. 

3. The number of certain types of operations such as the 
take-offs of long stage length aircraft may be limited. 
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A. Cerunlu runways lay only bo used during certain 

periods or a liil cation on the number of operations 
per period for any given runway may bo specified* 

5. Certain trajectories that correspond to specified 
ground tracks may only be used during certain periods 
or a limitation on the number of operations per period 
for any given ground track may bv specified. 

6. The total number of operations per time period may be 
constrained. 

Those flight limitations are now formulated mathematically: 

1, The annoyance as measured by some community response index 
for a given area, A, must not exceed some critical value, 

RI, < V, A = 1 NA (A. 2) 

For example if NEF -is selected as the connnunity response index, 
then (3.3) would appear as a function of the decision variables 

US 

NI NK 

10 log ^ T, £ 

i°l k=l 

-8.8) + EPNL^j^ + 10 logj^o^^'ijl “ 88D/10) £ 

A = 1, ..., NA (A. 3) 

2. The number of aircraft type t allowed to operate in and out of a 
given airport is limited to N^ operations per day 


NJ 

+ 10 log + 16.67X^y2) 
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NI 


(A./O 


NK 

T. 

k^l 


NJ 

K 

j“l 


NL 

Z 



NJ 

T. 


ML 


tjii 


< t “ 1 


3. The number of take-offs with stage length s is limited to N, 

s 

operations per day 


NI NJ 

NL 

Z Z 

Z 

1=1 j=l 

£=1 


^ isja •- s 


The number of operations for runway r is limited to Nr per day 


NT NK 
T. T. K 
i=l k»l Jt'R 


NL 
K X. 


NI 

+ j: I 


NL 

I y, 


'.i 


jeR i°l 
1 


r « 1 , , , ♦ , NR 


0. 


The number of operations corresponding to ground track d is 
limited to N^ per day 


Z Z T. X 

i k £ 


NI 

+ 5’. 


NL 


V Y < \ 

ikd£ , T id£ -» d 

a=l £=1 


d « 1 , . , . , N J 


6. The total number of operations in time period p must not exceed 
N 


NI NK NvT 


NI NJ 


Z Z Z X., , + Z S y, . < N p=l, NI 


i=l k=l j=l 


i=l j=l 
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A. 2. 2 Aircraft Availability 

Only a limited number of the various types of aircraft servicing an 
airport will be available during each time period for either landiiig or 
departure. This may be expressed analytically as 
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(A. 5) 


(A. 6) 


(A. 7) 


(A.8) 




i “ 1. . 
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(A. 9) 
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k“l j^l 


< NX 
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i » 1, ..., NX; S. “ 1, ..., NL 


(A. 10) 


A. 2. 3 PasBcnger - Aircraft Demand 

Passenger demand may be established for any given airport. Then 
passenger demand may be translated into aircraft demaiid. Such demands may 
take the form of requiring at least Nj.gjp type t aircraft operations of 
stage length s, utilizing ground track d during time period p. 


X 


tsdp 


> N 


tsdp 


t = 1, . . . , NX; s ° 1, . . NK; 
d = 1 , . . . , NJ ; p “ 1 , . . . , NL 


(A. 11) 


This is one of the simplest ways in which passenger demand may be accounted 
for. More elaborate demand relationships could be derived and utilized if 
desired. 

Collectively Equations (A. 2) through (A.H) define a mathematical model 
for this research. This model may be classified as a nonlinear Integer 
mathematical programming model. Solution techniques for such a model are 
discussed in Section IV and Appendix B. 


APPENDIX B 


Original Solution Algorithm 


Reference (7) 


Appendix A for the obj ^tive function form need, 


Step 1 


NA 


minimine « Z (P./^HIOKS. 

A-1 ^ ^ 
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(B,l) 


subject to all linear constraints where, 

EPNL 


N’l NK NJ ( 
S ° E E T. {10 
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letting, 
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and d^j. 
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This Av’ill provide a feasible solution to the problem but in no sense 
guarantees an optimum. Instead of the function given in Equation (4.1), 
the objective function should be 


N 

minimize 2^ = S (P^/P) (10) (log q(S^) - 1) 
A=1 


(B.4) 
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liov;cver the ourrognte objective func on has provided very good nolutions 
to the origiocl problem for the few ■ ‘sample problemo solved. 

The reason for such a ovirrogato objective function is that it is linear 
in the decision variables and subject to linear constraints, hence corresponds 
to a linear programming problem. With sophisticated computer implementation, 
linear programming solution techniques are capable of solving very large 
problems. For example, problems involving several thousand linear constraint 
equations and tens of thousands of variables arc within the realm of possibility 
for the very efficient linear programming solution algorithms Implemented on 
modern computers. However it should also be pointed out that to produce 
solutions to such large problems requires extensive efforts in data preparation, 
manipulation of the linear programming computer code, and interpretation of 
computational results . 

Step 2 

Obtain a truncated Taylor series expansion about tlvo solution point 
from Step 1, say _S* 

u(S) “ + VZ^(S*) (S “ S*) (B.5) 

Now minimize (4.5) subject to the original linear constraints. This 

corresponds to solving another liiaear programming problem. Denote the 

1 

solution to this linear program (4.5) as £ . Since u(^) is constructed from 
N 

the gradient of Z at an Improved solution point can be secured only 

if u(£^) < u(^*). This will not guarantee that Z^(£^) < Z^(£*) unless _S^ 

1 

is in the immediate neighborhood of However, given u(£ ) < u(^*) there 

2 1 
must exist a point, say £ , on the line segment between and £ such 

that Z^(S“) < Z^(S^’0. To determine ^ one solves 


rainimiao + a(S^ »- _S* 

a =. 


+ y(|.^ - 


S-'O) “ minlralac Z^(S* + a(S^ >* 

0<a<l ~ “* 


(B 


2 

Set S£ equal to S and repeat Step 2 an many tiroes as required to obtain 
the stopping condition, u(S_^) ^ u(£*) . At this point no further iroprovement 
is possible. 
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ALAG 1 


1, PURPOSE: To minimize a function F (x) “ f (X^, X^) 

subject to both equality and inequality constraints. 
Derivatives of all functions must be supplied in .a 
user subroutine entitled ALAGB (see item 5) . An 
Initial estimate of the solution (not necessarily 
feasible) must be specified. This computer program 
is developed from algorithm of section 

2. USE: CALL ALAGl (N,M,K,X,EPS, AKMIN, DFN, MAXFN, IPRl, 

1PR2, IW, MODE) 

N An INTEGER set to the number of variables 
n (N ^ 2) . 

M An INTEGER set to the total number of 

constraints m (M ^ 1) , 

K An INTEGER set to the total number of 

equality constraints k. 

X A REAL array of N elements in which the initial 

estimate of the solution must be set. ALAGl 
returns the solution x in X. 

EPS A REAL array of N elements, in which the 

tolerances for the unconstrained minimizations 
must be set. EPS (I) should be set so that 
EPS (I)/X (I) = AKMIN, roughly speaking. 
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AKMIN 


A REAL number in which the relative error 
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tolerance required in the constraint residuals 
must be set. ALAGI will exist when max{ | c.j, (x) | / 
scaling factor for ± AKMIN for the active 
constraints {i}. 

DFN A REAL number in which the likely reduction in 
F(x) must be set. This is done in the same 
way as for QN^-JTA - see the QNWTA description. 

MAXFN An INTEGER in which the maximum number of calls 
of ALAGB on any one unconstrained minimization 
must be set. 

IPRI An INTEGER controlling the frequency of printing 

from ALAGI. Printing occurs every IPRI iterations, 
except for details of increases to the c^ which 
are always printed. No printing at all occurs 
(except for error diagnostics) if IPRI = 0. 

IPR2 An INTEGER controlling the frequency of printing 
from QMTA. IPR2 should be set as described in 
the QNWTA documentation. 

IW An INTEGER giving the amount of storage available 
on C0MM0N/ALAGL/W(.) . Set to 2500 unless wishing 
to change the restrictions (see Section 5) . 

MODE An INTEGER controlling the mode of operation of 
ALAGI. If any positive definite estimate if 
available of the hessian matrix of the penalty 
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function, sa> jMODEj ° 2 or 3, otherwise set 
|M0 De 1 = 1 (see QNOTA description). If 
estimates of the and 0^ parameters are 
available (see item 8) set MODE < , otherwise 
set MODE >0. A normal setting for a one-off 
job with no information available is MODE “ 1. 

3. LABELED COMMON AREAS: 

Certain labeled COMMON areas must be declared and set on entry to ALAGl. 

C0MM0N/ALGAGE/C(150) Set scale factors (>0) for the constraints in 

C(l), C(2) ,G(M). Choose the magnitude of 

these scale factors to give an indication of 
the constraints evaluated about the initial ap- 
proximation X, If any constraints are violated 
by an amount greater in modulus than that which 
is set, then the setting is increased accordingly. 
These scale factors are transferred to C(M+1), 


C(M+2) C({aM) by ALAGl. 

C0MM0N/ALAGF/GC(25, 50) Set the derivatives of any linear constraints on 

entry rather than in ALAGB. This is the most 
efficient and the numbers are not disturbed. The 
manner of setting is described in item 4. 

If MODE < 0 is used, then set the parameters 
0^ in T(l), T(2) , . . . ,T(M) and the 
parameters ’’’’^m T(M+1) ,T(M+2) , . . , ,T(2M) 

The meaning of these parameters may be found in 
section of this report. 


C0MMON/ALAGG/T(150) 
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C0MM0N/ALAGI/G2P (325) 


If |mODE| “ .! or 3 set the estimated hessian 
matrix of the penalty function in G2P(1),..., 
G2P(N‘ (N+l)/2) . The manner of setting is that 
described in QNWTA under the heading MODE. 

Local storage for ALAGl is through labeled COMMON areas. These 
have been set on the assumption that N £ 25 and M £ 50. If it is desired 
to remove either or both of these restrictions, then it is necessary to 
increase the storage available in some or all of these areas. This can 
be done by defining the named COMMON areas in the users MAIN with the 
increased storage settings, in which case the extra storage will be 
effective throughout the whole program. The complete list of labeled 


COMMON used by ALAGl and the corresponding values of N and M are as follows. 


COMMON/ALAGC/F , M, K , IS , MK , NU 
" D/G(50) 

" E/C (150) 

" F/GC(25,50) 

" G/T(150) 

" H/GP(50) 

" I(G2P(325) 

" J/V(50) 

" K/W(150) 

" L/W(2500) 

" M/ZZ(100) 

" N/LT(100) 


independent of N and M 
2N 
3M 
N,M 
3M 

y (y = max(M,N)) 
N- (N+D/2 
y 

3y 

2 

y 

2y 

2M 


4. ACCUMCY: Tills iterative alg-Tithm terminates normally when 

the following convergence condition is met; 

nuiK { jc^(x) I /scaling factor for c^} £ AKMIN for 

i an element of the set of active constraint indices, 

A diagnostic message for abnormal termination is 
printed when the program is unable to achieve the 
requested accuracy. This may be due to (i) a 
mistake in programming ALAGB, (ii) there is no 
feasible point (in which case ->• ® and c^ ->• constant 
^ 0), (iii) EPS has been set too large relative to 
AKMIN, (iv) the problem is too ill-conditioned. 

OTHER ROUTINES: ALAGl requires the use of ALAGB, ALAGZ, BQDMA, 

MULDA, MULDB, MULDE, and QNWTA 

ALAGB: USER SUBROUTINE The user must define a subroutine headed by 

SUBROUTINE AUGB(N,M,X) 

REAL X(l) 

COMMON/ALAGC/F 
COMMON/ ALAGD/G (50) 

COMMON/ ALAGE/C (150) 

GOMMON/ALAGF/GC (25 , 50) 

This subroutine takes the vector X and sets 

(1) F(x) in F; (2) c (x) , . . . , c (x) in G(l) , . . . ,G(M) ; 

X m 

(3) (DF/-„ ) |x in G(1),...,G(N); 

ax^,...,3F/g^^ 

n 

(4) Oc^/gj. ... ac / )F in GG(1,I),... 

1 ’ ' * ' ’ i 9Xn 

GG (N,I) for I = 1,...,M. 
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ALAGii: This subroutine evaluates tlie augmented function comprised 

of the original objective function and penalty terms that is to be 
optimized. 

SUBROUTINE ALAGZ (N, X, PHI, GPHI) 

N and X as previously defined, 

PHI is the value of the augmented function evaluated at X. 

GPHI is the gradient of the augmented function evaluated at X. 

BQUMA: The purpose of BQDMA is to find the values that minimize a 

quadratic of n variables subject to upper and lower bounds on some or 
all of the variables subject to upper and lower bounds on some or 
all of the variables. 

The quadratic is defined by 

Q(X) = 1/2 X*^ AX - B*^X 
Subject to: 

1 \ 1 2U^ i = 1,.. .N. 

'■UBROUTINE BQDMA (N,A,IA,B,BL,BU,X,Q,LT,K,G) 

N an INTEGER which must be set by the user to the number 

of variables. 

A a REAL, two dimensional array, each dimension at least 

N; the elements in the upper triangle A(I,J) l£J^N must 
be set by the user to the corresponding A^^ in (1) , and 
will remain untouched by the subroutine. Elements 
A(I,J) N^I>J are used as working space. 
lA ..u INTEGER giving the first dimension of A in the 

statement which assigns space to A. 

B a REAL array of at least N elements. The user must set 

b(I), B is not oven^ritten by BQDMA. 


ALAGZ: This? subroutine evaluates t!io augmented function comprised of 

the original objective function and penalty terms that is to bo optimized. 
SUBROUTINE ALAGZ (N, X, PHI, GPHI) 

N and X as previously defined. 

PHI is the value of the augmented function evaluated at X. 

GPHI is the gradient of the augmented function evaluated at X. 

BQDMA: The purpose of BQD'MA is to find the values that minimize a 

quadratic of n variables subject to upper and lower bounds on some or 
all of the variables. 

The quadratic is defined by 

Q(X) = 1/2 X^ AX - B*^X 
Subject to: 


BL, < X, < BU. 
1 i ™ 1 


i = 1, , . . ,N. 


SUBROUTINE BQDMA (N,A, IA,B,BL,BU,X,Q,LT,K,G) 

N an INTEGER which must be set by the user to the number 

of variables. 

A a REAL, two dimensional array, each dimension at least 

N: the elements in the upper triangle A(I,J) l£JfN must 
be set by the user to the corresponding A^j in (1) , and 
will remain untouched by the subroutine. Elements 
A(I,J) N^I>J are used as working space. 
lA an INTEGER giving the first dimension of A in the 

statement which assigns space to A. 

B a REAL array of at least N elements. The user must set 

B(I). B is not overwritten by BQDMA. 


a REAL array pf at 1> ist N elements. The user must 

th 

set LL(I) to the lower bound on the *■ variable. 

If the bound is non-existent, set it to a very small 
number like -1E75. BL is not ovenn'ilrten by BQDMA. 
a REAL array of at least N eJements. The user must 
set BU(I) to the upper bound on the variable. 

If the bound is non-existent, set ti to a very large 
number. BU is not overwritten by BQDMA. 
a REAL array of at least N elements. BQDMA returns 
the solution in X(I) . 

a REAL variable in which BQDMA returns the solution 
value of the quadratic. 

an INTEGER array of at least N elements, set by 
BQDMA to a permutation of the integers 1,2,...,N 
(see K and G below) 

an INTEGER set by BQDMA to the number of free 
variables at the solution (those not on their bounds). 
These are the variables LT(1) , LT(2) , . . . ,LT(K) . 
a REAL array of at least 3*N elements. G(l) , . . . . ,G(N) 
are set by BQDMA to the gradient evaluated at the 
solution point. G is indirectly addressed so that 
G(l) contains the gradient with respect ot the LT(I) 
variable, wlience G(l) , . . . ,G(K) will be found to be 
zero. G(N +1) , . . . ,G(3*N) are used by BQDMA as working 
space. 
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MULDA is a subroutine for use iti problems which involve the 

T 

addition or subtraction of rank one matrices o to positive 

definite or semi-definite symmetric matrices A stored in factored 
T 

form A == LDL , such that the resulting N xN matrix 

* A ^ T 

A « A + g zz 

is also kno\'m to be positive definite or serai-definite. Note that L is 
lower triangular with “1, and D is diagonal with d^ ^ 0, 

SUBROUTINK MULDA (A, . . L', IK, MK , EPS) 

A A REAL one dimensional array of N*(N+l)/2 elements 

T 

in which the matrix A=LDL must be given in factored 

form. The order in which eleravints of L and D are 

s tored xs d- > i * ^‘x i * * * * > o s 

X 4X NX 4 , • . • , ^^2 ^ * * * * ’ * 

dx, ij^x, M TjdvT* I'he. factors of the matrix 
N-1 N,N“1 M 

A “ A + o zz will overwrite those of A on exit. 

N An INTEGER (N^^l) which must be set to the dimension 

of the problem. 

Z A REAL one dimensional array of N elements in which 

the vector z must be set. The array Z is oven^ritten 
by the routine. 

SIG A REAL variable in which the scalar a must be set. 

SIG is not restricted to +. , but if SIG<0 then it 
must be known from other consideratii'.., ;.at A is 
positive definite or semi-definite, apart from the 
effects of round-off error. 
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A REAL array of N eli neats. If SIG>0 then W is 
not used, and the name of any one dimensional array 
can be inserted in the calling sequence. If 
SIG>0 then W is used as work space. In addition for 
SIG^O it may be possible to save time by setting in 
W the vector v defined by Lv^a. The ways in which 
this can occur are described under MK below. 

An INTEGER to be set that |ir 1 is the rank of A. 

If the rank of A is expoctud to be different from that 
of A, set IRj$,0. On exit from MULDA, IR(^O) will 
contain the rank of A. 

An INTEGER to be set only when SIG<0, as follov'js. 

If the vector v defiiied by Lv!=z has not been calculated 
previously, set If MULDA has been used previously 

to calculate A ^ z, then v is a by-product of this 
calculation and is stored in the W parameter of MULDE. 

In this case transfer v to the W parameter of MULDA 
and set MK“1. If z has been calculated as z = Au for 
some arbitrary vector u using MULDD, then again v is 
a by-product of the calculation and is available in the 
W parameter of MULDD. In this case (or any other in 
which V is knom) set v in the W parameter of MULDA 
and set MK=2. 

A REAL variable to bo set only when SIG<0 and A .is 
expected to have the same rank as A. In certain ill- 


conditioned cases a non-zero diagonal element 

o£ D might become so small as to bo indeterminate, 

IVo courses of action are possible. One is to 

introduce a small perturbation in order that A 

keeps the same r*ank as A. This is the normal course 

of action and is achieved by setting KPS equal to 

the relative machine precision e. The other course 

** 

of action is to let the rank of A be one loss than 
the I'ank of A, This is achieved by setting EPS 
equal zero. 

MULDB - factorizes a positive definite syramotrle matrix given in 
A. This matrix is then used in MULDA. 


SUBROUTINE MULDB (A, N, IR) 
A 


N 

IR 


Must contain the elements of A in the order 

’‘^21’ ■ • ’ ’**N1’‘‘*22’**32’ • ‘ •‘'*N2’ ' ‘ ’ ’*Vl,N-l’**NN5 
that is as successive columns of its lower triangle) . 
One exit A will be overwritten by the factors L and D 
in the form described in MULDA. 

Order of the matrix A. 

An INTEGER set by MULDB to the rank of the factori- 
zation. If the factorization has been performed 
successfully IR=N will be set. If IR<N then the 
factorization has failed because A is not positive 
definite (possibly due to round-off error). In this 
case the factors of a positive semi-definite matrix 
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of rank IR will be f'?und in A. However Che 
rcoulta of thio calculation are unpredictable, 
and MULDD should not be used in an attempt to 
factorize positive serai-definite matrices. 

MUhDE calculates the vector z* » z where A is in factored form 
SUBROUTINE >RJhUE (A, N, E, W, IR) 

A Must be set in factored form. 

N Order of the matrix A. 

7 , A REAL array of N elements Co be set to the vector z. 

On exit Z contains the vector z* “ A ^ z. 

W A REAL array of N elements which is set by MULDE 

to be vector v defined by Lv°z, If this vector is 
not of interest, replace W by Z in the calling 
sequence to obviate the need Co supply extra storage. 

IR An INTEGER which must be set to the rank of A. 

QNTWA finds the minimum of a function F(x)of several variables 
given Chat the gradient vector can be calculated. This routine is based 
upon a quasi-Newton method described by Fletcher in (F8) . 

SUBROUTINE QNWTA (FUNCT, N, X, F, G, H, W, DFN, EPS, MODE, MAXFN, 

IPRINT, lEXIT) 

FUNCT An IDENTIFIER of the users subroutine. 

N An Integer to be set to the number of variables (N ^ 2) . 

X A REAL ARRAY of N elements in which the current estimate 

of solution is stored. An initial approximation 
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rauot be set in X on i>ntry to QNWTA and the bent 
eatiraate obtained will be returned on exit. 

A REAL number in which the best value of F(x) 
correaponding to X above will be returned. 

A REAL ARRAY of N elements in which the gradient 
vector corri“;p*n; : ng to X above will be returned, 
f it) be set on entry. 

A REAL ARRAY of N*(N+l)/2 elements in which an 
estimate of the hessian matrix is stored. The 
mati’ix is represented in the product form LDL 
where L is a lower triangular matrix with unit 
diagonals and D is a diagonal matrix. The lower 
triangle of L is stored by columns in H excepting 
that the unit diagonal elements are replaced by 
the corresponding elements of D. The sotting of 
II on entry is controlled by the parameter MODE. 

A REAL ARRAY of SJ'fN elements used as working space. 

A REAL number which must be set so as to give QNI^TA 
an estimate of the likely reduction to be obtained in 
F (x) . DFN is used only on the first iteration so 
an order of magnitude estimate will suffice. The 
information can be provided in different ways 
depending upon the sign of DFN which should be set 
in one of the following ways: 


DFN>0 


EPS 


MODE 


the 8G! ting of DFN itself will be 
taken as the likely reduction to be 
obtained in F(x). 

DFN“0 it will be assumed that an estimate of 
the minimum value of F(x) has been set 
in argument F, and the likely reduction 
in F(x) will be computed according to 
the initial function value. 

DFN<0 a multiple |dfk| of the modulus of the 
initial function value will be taken as 
an estimate of the likely reduction. 

A REi\L ARRAY of N elements to be set on entry to 
the accuracy required in each element of X. 

An INTEGER which controls the setting of the initial 
estimate of the hessian matrix in the parameter H. 

The following settings of MODE are permitted. 

M0DE'=1 An estimate corresponding to a unit 
matrix is set in H by QNMTA. 

M0DE==2 QNWTA assumes that the hessian matrix 

itself has been set in H by columns of 

its lower triangle, and the conversion 
T 

to LDL form is carried out by QNWTA. 

The hessian matrix roust be positive definite. 

M0DE=3 QNWTA assumes that the hessian matrix has 
been set in H in product form. This is 
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convcnilint when using the H matrix 
from .one problem as an initial 
estimate for another, in which case 
the contents of H arc passed on 
unchanged . 

MAXFN An INTEGER set to the maximum number of calls of 
FUNCT permitted. 

ItlPUINT An INTEGER controlling printing. Printing occurs 

every [iPRINT] iterations and also on exit, in the form 
Iteration No, No of calls of FUNCT, lEXIT (on 
exit only) . 

Function value 

X(1),X(2),...,X(N) 8 to a line 

G(I) ,G(2) , . . . ,G(N) 8 to a line 

The values of X and G can be suppressed on inter- 
mediate iteratior;' by setting IPRINT<0. All 
intermediate printing can be suppressed by setting 
IPRINT-MAXFN+1. All printing can be suppressed by 
setting IPRINT=0. 

lEXIT An INTEGER giving the reason for exit from QinWA. 

This will be set by QNWTA as follows: 

IEXIT=0 (M0DE=2 only). The estimate of the 

hessian matrix is not positive definite. 
IEXIl-1 The normal exit in which |dX( 1) [ <EPS (I) 
for all 1=1, 2,... N, where DX(I) is the 
change in X on an iteration. 
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IEXIT°2 G^DX; U. Not puf,sl.ble without rounding 
error. Probable cause is that EPS is 
set too small for computer word length. 


IEXIT=3 FUNCT called IMAXFN times. 



APPENDIX D 


COM’UTATIONAL RESULTS POR ALAGl ALGORITIM 







COMPUTATIONAL RESULTS FOR ALAGl OPT' MIEATION ALGORITHM 


Nil a 0.15276 


Numb 01 - of 
Opera fions 


Aircrof C 

Arrival Departur e 727 I 737 


Stape Track 
I.enpch Number 


Period 

Night 


APPENDIX E 


SanHitivicy Analyses P’or Example Airport 
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The mathemaclcal model employed *or nearly all the sensitivity analyses 
is presented below. The only result! that were not derived through the use 
of this model are those in Tables ElO and Ell. These results are for a 
"minimax" objective function described in section four. 

The mathematical model for this second example problem is 

31 

minimize K (P^/P) (10) (log^gCS^) - 1) (E.l) 


subject to: 


NI 

^^111 ^1121 \ l 31 ^ 


> 3 


(E.2) 


Nl 

^^211 '\ 231 ^ 


> 3 


(E.3) 


NI 

*^^1311 ^321 ^ 331 ^ 


> 2 


(E.4) 


NI 


> 6 


(E.5) 


NI 

/' ^^1241 ^251 \261 ^ 271 ^ -- 

X^X 


(E.6) 


NI 


‘‘^1341 ■** ^351 ^361 ^371) ^ 5 


(E.7) 


NI 

^^1131 -*-w 

x=l 


> 2 


(E.8) 


NI 

^^1281 ^ 1291 ^ 


> 4 


(E.9) 
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V 


NI 

^^1381 ^'1391^ 


> 3 


NI 

i“l ^^2A2 ■** ^252 ^1262 '’‘ ^1272^ ~ ^ 


NI 

^^1182 \l92^ 


> X 


NI 

V (Y + Y + Y I 
/' ^ ill 131-’ 

i“l 


> 20 


NI 

^'’^141 '"'isi \ei 


> 14 


NI 

^\81 ^i91> 


> 3 


NI 

F (Y + Y' + Y I 
^’‘112 " l 22 ^ 32 ^ 


> 8 


NI 

V (Y + Y ■) 
^^182 ’^192^ 


> 1 


NJ 

S Ijl 


j=l 


< 11 


NJ 
Z Y 
j=l 


2jl 


< 28 


Nj 
^ Y 

j=l lj2 


< 3 


- 60 - 


-dk^ -dfcfe 4 j 


(E.IO) 


(E.ll) 


(E.12) 


(E.13) 


(E.14) 


(E.15) 


(E.16) 


(E.17) 

(E.18) 


(E.19) 


(E.20) 


CR.21) 


N.T 

K 

j-1 



< 7 


NK 

V 

k«l 


NJ 

T. 


^2kjl 


< 9 


NK NJ 

V V 

k«l j«l 


^'2kjl 


< 32 


NK NJ 

V V 

k=l 


^lkj2 


0 


NK 

j; 

k=l 




< 1 


NI NI 

}: Y }: 

i-l ^ i=l 


NK 

V 

k=l 


^ikj 4 il 


1 , • • • ) NL; 
1, .... NJ 


Table El details the sensitivity analyses conducted for the 
airport whose mathematical model was just presented, 
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(E.22) 


(E.23) 


(E.24) 


VK.25) 


(E.26) 


application 




TAB! El 


Demand Ve?iations 


Right Hand 

Side of Constraint 

Case 1 

Case 2 

Case 3 

Case 4 

E2 

5* 

3 

2* 

3 

E3 

5 * 

3 

2* 

3 

E4 

3* 

2 

1* 

2 

E5 

10*v 

6 

4* 

6 

E6 

16* 

11 

9* 

12* 

E7 

8* 

5 

4* 

5 

E8 

3* 

2 

1* 

3* 

E9 

6* 

4 

3* 

4 

ElO 

5* 

3 

2* 

3 

Ell 

1 

4* 

1 

0* 

El 2 

1 

4* 

1 

0* 

E13 

36* 

20 

14* 

28* 

E14 

24* 

14 

9A 

14 

E15 

5* 

3 

2* 

4* 

E16 

8 

9* 

5* 

0* 

E17 

1 

4* 

1 

0* 

E18 

20* 

11 

9* 

14* 

EL9 

48* 

28 

21* 

35* 

E20 

3 

5* 

2* 

0* 

E21 

7 

10* 

5* 

0* 

E22 

18* 

9 

7* 

10* 

E23 

44* 

32 

24* 

33* 

E24 

2 

6* 

2 

2 

E25 

1 

4* 

1 

1 


Designate changed values. 



TABLE *:2 


W 


n 


COMPUTATIONAL RESULTS FOR DEMAND INC 'ASE FOR DAY FLIGHlb 


Nil a 0.13026 
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The various cases correspond to the ' JlXowing variations. 


Case 1 - Demand increase for day flights 
Case 2 “ " " "night " 

Case 3 ~ Uniform decrease in demand 

Case -fi - All flights are conducted during the day 

The computational results are compiled in Tables E2 through E5. 

Another type of sensitivity analysis conducted for the given application 
was variations in the type of aircraft employed for all operations. The 
specific cases considered were: 

Cate 5 - all operations were by 737 aircraft 

with SAM 
with SAM 

Tables E6 through E9 provide the corresponding computational results. 

Table KIO describes the results obtained when the objective function 
was changed from the one originally given in aquation (El) to that of mini- 
mizing the maximum (minimax) annoyance for the various population segments 
considered. The use of such an objective formulation should result in a 
dispersion of annoyance more uniformly to the population areas surrounding 
the airport. The solution of a population weighted NEF minimax objective 
function is given in Table Ell, 


Case 6 - " " " " " 

Case 7 - " " " " 727 " . 

Case 8 - " " " " " " 





















































































TABL E4 
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XA ',K E6 

COMPUTATIONAL RESUI.TS FOR ALL 737 AI CRAFT 









































































































































































TABL E9 


COMPUTATIONAL RESULTS FOR ALL 727/SA^' AIRCMFT 


Nil ° 0.17922 


Numb or of 
Oiuiracloi'iH 

Arrival Departure 

Aircr 

727 

aft 

737 

Stage 

Length 

Track 

Number 

Time 

Day 

Period 

Night 

ZO 



X 



t 

X 


6 

K 


X 



1 


X 

1^ 



X 



4 

X 


3 

K 


X 



9 



1 

X 


X 



9 


X 

3 


,X 

X 


1 

f 

X 


3 

— 

X 

X 


2 

\ 

X 


t 


X 

X 


3 

\ 

X 




v> 

/\ 

X 


1 

A 

1 

✓\ 


11 


X 

X 


1 

4 

X 


s 


X 

X 


3 

4 

X 


1 

r 

X 

X 


2 

4 


X 

2 


X 

X 


\ 

9 

X 


4 


X 



2 

3 

X 


3 


X 

m 


3 

3 

X 


1 


X 

X 


1 

3 


X 




m 


















im 









■ 
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■ 
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TAB ElO 

COMPUTATIONAL RESULTS FOR MINIMIZE M.' TIMUM ANNOYANCE 


Nil = 0.15968 


Number of 
Operations 

Arrival Dcviarture 

Aircx 

727 

•aft 

737 

Stage 

Length 

Track 

Number 

Time 

Day 

Period 

Night 

9 

X 1 



X 


2 

X 


11 

X 


X 



3 

X 


3 

X 


X 



3 

• 

X 

1 

(> 

K 



X 


3 


X 

1 

X 



X 


C> 

X 


13 

X 

, 


X 


7 

X 


3 

X 



X 


B 

X 


1 

X 



X 


3 


X 

3 


X 

X 


t 

1 

X 


3 


X 

X 


2 

1 

X 


1 


X 


X 

3 

3 

X 


1 


X 

X 


- 1 

4 

X 


S 


X 


X 

1 

4 

X 


11 


X 


X 

2 

4 

X 


. S 


X 


X 

3 

4 

X 


1 


X 

X 


Z 

4 


X • 



X 


X 

I 

3 

X 




X . 


X 

Z 

3 

X 

* 

3 - 


. 


> X 

5 

3 

X 


1 


X 


X 

t 

9 


r 

X 

i” " 









k 









' * 
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OriRinal Objective F -netion Mininuix 

Nil Contribution NEF Nil Contribution 


Objective Function 


.0A9 
.176 
.137 
.438 
.738 
1.038 
.143 
.032 
.019 
.705 
.581 
.161 
. 300 
. 362 
.352 
.378 
.926 
.338 
.140 

.542 

1.140 

.867 

1.332 

,416 


8.57 

11.41 

10,07 

15.82 

23.99 

27.05 

20.66 

12.43 

4.87 

19.36 

20.90 
20.14 

31.66 
23.75 
17.16 
18.50 
18.98 
23.58 
39.19 

24.40 
23 . 62 

26.67 

35.91 

28.40 


.295 
.676 
.567 
.896 
1.071 
.498 
.084 
.061 
1.062 
1 . 004 
.263 
.374 
.147 
.670 
.697 
1.723 
.527 
.172 

.823 

1.371 

.520 

.836 

.209 


14.86 
21.16 
17,63 
25.61 
27.33 
31.36 
19 . 25 
12.20 
22.48 
25.28 

23 . 98 

33.98 
16.94 

21.91 

23.14 
2 U76 

. . * ■» i 

41.91 

28 . 04 

25.14 
22.42 
31 . 18 
22.54 


.277 
.169 
.167 
1 . 054 


14 . 24 
12.92 
22.04 
33.79 


19.59 

18.52 

19.49 

22.41 


12.979x10 


15.968 X 10‘ 

















































































































